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The topological properties of a multiband topological superconductor in two dimensions are stud¬ 
ied, when the latter is obtained by introducing electron pairing in an otherwise topological insulator. 
The type of pairing, doping and Fermi surface topology play an essential role. Considering the An¬ 
dreev reflection problem, we use a previously developed quantum waveguide theory for multiorbital 
systems and find that, when the Fermi surface has several pockets, this theory retrieves the correct 
number of Majorana fermion states as predicted by the topological index. By varying band structure 
parameters, the Fermi surface topology of the normal phase can be made to change, whereby the 
number of Majorana modes also varies. We calculate the effect of such transitions on the Andreev 
differential conductance. 

PACS numbers: 71.10Fd, 71.10. Pm, 71.10. Li, 73.20.-r, 74.45.+C, 74.20.Rp 


I. INTRODUCTION 

Recent interest in non-trivial topological properties of 
insulator^ has spurred intensive research on band mod¬ 
els displaying non-trivial topologjffl. In the case of super¬ 
conductivity, non-trivial topology is associated with the 
emergence of zero energy excitations which are their own 
antiparticles, or Majorana fermions (MFs). Condensed 
matter therefore provides the ground for the realization 
of such long sought exotic particles 3 4 . 

Early theoretical models of two-dimensional topologi¬ 
cal superconductivity consider p + ip pairing in an oth¬ 
erwise trivial band. In such a framework, the emer¬ 
gence of MFs in spinless fermion models was shown to 
be related to topolog)®. Non-abelian MF’s were shown 
to arise at vortex cores in a model of spinfull fermions 
with spin triplet p + ip pairing and where the two spin 
components effectively decouple®. It is this concept of 
MF that we shall deal with, in this work. Besides their 
interest from the point of view of fundamental physics, 
the MF’s non-Abelian exchange statistics is their most 
remarkable property and suitable for quantum compu¬ 
tation applications 6 . It was early realized that an ef¬ 
fective p + ip superconductor arises from single band 
Dirac electrons complemented with s-wave pairing® A 
promising venue for the construction of topological su¬ 
perconductors is the proximity coupling of a topologi¬ 
cal insulator to a superconductor^®. The possibility of 
engeneering topological superconductors using the sur¬ 
face of a three-dimensional TI or a two-dimensional semi¬ 
conductor in proximity to a s-wave superconductor was 
discussed 8 but multiorbital effects were not considered. 
Non-conventional superconductivity and Andreev reflec¬ 
tion were also studied®. We note that topological ma¬ 
terials are necessarily multiorbital. The consideration of 
p + ip (or other) pairing in a single band model is there¬ 
fore a simplification. 

More recently, there has been interest in multior¬ 


bital systems. The orbital degree of freedom (pseu¬ 
dospin) is present in a lattice symmetry transforma¬ 
tion, namely, that of lattice inversion, and the char¬ 
acterization of pairing term as having even/odd parity 
under inversion. Such investigations have been moti¬ 
vated by the three-dimensional topological superconduc¬ 
tor CuzB^Sea, which has two orbitals per lattice cell. 
Some theorems on the topological indices to be expected 
for various superconductor models, by taking their sym¬ 
metries into account, have recently been established. In 
many cases, the topology of the Fermi surface itself is 
important. For instance, under the assumption of lat¬ 
tice inversion symmetry, time-reversal invariance (TRI) 
and odd parity pairing, three-dimensional superconduc¬ 
tors are topological if they possess an odd number of 
Fermi surface pocketf 1112 . Models of two-dimensional 
superconductors with a pseudospin degree of freedom 
have been proposed recently, concentrating on the case 
of nodeless odd parity pairing in TRI superconductors. 
In this case non-trivial topology requires spin -orbit cou¬ 
plings non-diagonal in the pseudospin channel- 13 ^ 14 ! 

we introduce below a model for a superconductor 
which has two orbitals per lattice site. For the case of 
diagonal (in pseudospin space) pairing, we find that the 
topological properties of the kinetic energy have no influ¬ 
ence on those of the superconductor. The latter are deter¬ 
mined, instead, by the pairing symmetry and the Fermi 
surface (FS) topology. We also study the Andreev prob¬ 
lem at the normal/superconductor (N/S) boundary, as¬ 
suming the normal metal to be single band, using the pre¬ 
viously established quantum waveguide theory (QWT). 
QWT was originally developed to address the And reev 
reflection problem in heavy-fermion superconductors 15 16 
and, later on, iron pnictide superconductors 1 7 3S5I. It can 
be applied when the N/S interface is parallel to one 
the primitive vectors of the superconductor and normal 
metal. For other interface angles, other approaches have 
to be employed. QWT also ignores the detailed micro- 
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scopic behavior of hopping parameters at the interface. 
Refinements of this theory which consider fully micro¬ 
scopic behavior of hopping parameters can be found in 
Ref sPSI. 

While the Blonder-Tinkham-Kaplwjck (BTK) 
theory 23 predicts a zero energy Andreev bound state 
(ABS) for each FS pocket, QWT accounts for the quan¬ 
tum interference effects between different FS pockets. As 
we shall see below, the ABS predicted by QWT correctly 
reconciles the Andreev problem with the topological 
index. 


II. MULTIORBITAL SUPERCONDUCTOR 
A. Pairing in a two-dimensional TI 


this case the eigenfunctions of Q can be written as a 
direct product 



Wave function normalization reads \u\ 2 + \v\ 2 = \a\ 2 + 
|/31 2 = 1. The amplitudes (a,/3) diagonalize the kinetic 
energy 

[h(k) ■ r + /io(k)ro] ^ ^ =£(k) ^ ^ , (6) 

and the BdG amplitudes ( u , v ) obey 

(7) 


The simplest version of superconducting pairing is that 
where electrons with opposite spins are related by time- 
reversal in the normal phase. The model includes at least 
two orbitals per lattice site (pseudo-spin) in addition to 
the spin degree of freedom. In the simplest version, the 
spin f electrons have kinetic energy S|(k) which is related 
to that of the spin f electrons by time-reversal, S|(k) = 
Sl(—k). If has an odd Chern number then the normal 
system serves as a model for a TI in two dimension^. In 
the case of two orbitals per lattice cell, we can write 

H t (k) =/i(k)-r + /i 0 (k)r 0 (1) 

where the Pauli matrices 7 y = i .2,3 act on orbital space 
(pseudo-spin) and r 0 denotes the identity matrix. We 
shall later use Pauli matrices cq, ti and r* operating on 
spin, particle-hole and Bogolubov-de Gennes amplitude 
(u v) spaces, respectively. 

The Bogolubov-de Gennes (BdG) matrix in the 
particle-hole basis (c|C|c|c|), where the field operators 
(:„ include the pseudo-spin degrees of freedom, takes the 
form: 



where 2 = diag(S|, S.Q. Keeping spin s z as a good 
quantum number, the pairing matrix A can have sin¬ 
glet ip(k)i<j 2 and triplet d~(k)eri components. The form 
of the the pairing matrix in pseudospin space remains to 
be chosen. We shall write 

A(k) = (A s il;(k)ia2 + A t d z (k)cr 1 ) tj , (3) 


Where we have defined the vector h' components as h! x — 
ih' y = A(k) and h' z = f(k). The two normal state bands 
are £(k) = ±|/r(k)| -f ho(k)T 0 . The Berry connection 
obtained from 0 is 

A = j GF '"A 1 k (w) 

= a' + a (8) 

where a' and a denote Berry connections associated with 
{u v) and (a /3), respectively. The total Chern number, 
C , is given by the line integral around the BZ, 


C = Yrf A ' Sk 


(j a' ■ Sk + (j a ■ dk 


(9) 

( 10 ) 


where the summatio n in ([ 9 ]) is over th e two negative BdG 
bands E± (k) = — y/^.(k) + |A(k)| 2 of the eigenproblem 

([ 7 ]). The summation in the second term of equation (10) 
is the Chern number of the two normal bands C(k) and 
therefore equates to zero. Hence we find that the topology 
of S| does not contribute to C because of cancelation 
in the sum over the normal state bands. It remains to 


analyse the first term in equation (10). Using Stokes 


theorem, the line integral of a' can be written as the flux 
of a monopole f2' = V x a' = h'/(2h 12 ), through the 
Bloch sphere: 


where the Pauli matrix Tj remains to be specified. Then 
H splits into two BdG matrices for (c|cj)) and (cj.c|) 
spaces. The former reads 

H = (/i(k) • t + /i 0 (k)r 0 ) r 3 

+ Re [A(k)] Tjri - Im [A(k)] Tjr 2 , (4) 

We next study the topological indices of 0- We first 
consider diagonal pairing in orbital space, Tj = tq. In 


c =^Y.I a - dS '*- (u) 

Here, h' = ^ReA(k), ImA(k), are the two 

monopoles’ curvature fields. The vector h! covers the 
Bloch sphere C times. The north (N) and south (S) hemi¬ 
spheres correspond to h' z = £ being either positive or neg¬ 
ative. If the chemical potential, which is included in /iq, 
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lies in the gap between the normal £± bands, then C = 0 
because h' stays always in one hemisphere. Nonzero C 
requires the chemical potential to intercept at least one of 
the normal bands (£_ or £ + ) and the N/S hemispheres are 
attained when h' z is outside or inside FS pockets. This 
links C to the topology of the FS. The other normal band 
does not contribute to C. The N and S poles are attained 
at BZ points where h' x ± ihf y = A^*^(k) vanishes. This 
motivates the choice of p + ip pairing below. A relevant 
example of a multiband sup ercon ductor believed to have 
p +ip symmetry is Sr 2 RuOj 24 1 25 1 

We recall that a theorem relating the topological in¬ 
dices of a superconductor to the FS topology has been 
established by Satd^for the case where the normal bands 
have inversion symmetry and pairing has odd parity. Our 
discussion above makes no requirement on inversion sym¬ 
metry or parity. In the case of time-reversal invariant 
single band spin triplet superconductors, the topolo gica l 
indices were also shown to be related to FS topology 2 ^. 

For other choices of Tj in equations |3])-([4]) the direct 
product form <[5j) is no longer valid. We shall address 
those cases below by computing winding numbers for TRI 
momenta where H is chiral. 


B. Model 


We write kinetic energy for f-spin electrons as in equa¬ 
tion 0 , where 

h x = sin k y , h y = — sin k x , 

h z = 2ti (cos k x + cos k y ) + 4t 2 cos k x cosk y , 

ho = — p — t\ {cosk x + cos k y ) . (12) 


We consider a spin triplet p + ip pairing: 


A(k) = d z (k)cri ® t 0 . 


(13) 


with d~(k) = A (sinfc x — isink y ). The pairing term (13) 
has odd parity under inversion. This symmetry property 

read^ r 1 A(-k)r 1 = -A(k). 

In the absence of Rashba or Dresselhaus spin-orbit cou¬ 
plings the Bogolubov-deGennes matrix in the particle- 

hole basis, ^(c|)(c|)(c|)(c|)^, splits into two 4x4 ma¬ 
trices. We consider now only the subspace ^(c-f)(c|)^. 


The parameter choice p = 0.6, t\ = 0.07, t 2 = —0.08, 
produces a FS with 3 pockets centered at (0,0), (0,7r), 
(■7 t, 0) in the Brillouin Zone (BZ), as figure [I] (left panel) 
shows. We calculate the topological index (Chern num¬ 
ber) following a method for multiband systemiPU and ob¬ 
tain C = +1, indicating that one MF exists. The energy 
spectrum for an infinite ribbon in the yy (or xx) direc¬ 
tion is shown in figure [2} where the MF is clearly seen at 
longitudinal momentum 7r. If we consider the Andreev 
problem for a N/S interface along yy, we expect the MF 
to be detected when the incident electrons have trans¬ 
verse momentum k y = n and thus traverse only one FS 
pocket, at (0,7r). 


By reversing the sign of the hopping parameter t 2 , a 
topologically trivial phase is obtained, with zero Chern 
number, and the FS now contains 4 pockets, as can be 
seen from the right panel of figure [l] and no MF’s should 
exist. 

One can also consider even parity p airi ng, which is 
realized by changing the r matrix in (13) to A(k) = 
c4(k)cri (gi T3. A direct computation of the ribbon spec¬ 
trum and Chern number yields the same topological in¬ 
dices as before. The ribbon spectrum remains gapped in 
the vicinity of k = 7r, displaying the same edge mode, 
similar to that in Figure [2] This case lies outside the 
conditions of the theorem proved above. 

From equations Q and ( |l3| it is clear that the Hamil¬ 
tonian for each subspace has the quiral properties 29 : 


riH(k)ri = -H( k) if k x = 0,7r, (14) 
r 2 .ff(k)r 2 = —H{ k) if fc y = 0,7r. (15) 


In either case one can rotate H to the off-diagonal form, 

U*H(k)U = ( ^ q ) . (16) 


The matrix U is composed of the column eigenvectors of 
rq or r 2 , for each of the chiralities (14)-(15), yielding 


Ai (k y ) — — a, + A, (17) 

A 2 (k x ) = — + ?A , (18) 


respectively. The phase of the determinant of A 1 ^ 2 ) accu¬ 
mulates an amount 2ttWu 2 ' ) , as its argument, k y r x \, goes 
from —7r to 7 r. The integer W is the winding number. 

For the choices t :] = ro ,3 in equation ([ 3 ]) and t 2 = 
—0.08 (topological case), we obtain W 2 (k y = 7r) = 1 and 
Wi(k x = 7r) = — 1. If t 2 = +0.08 (non-topological case), 
then all Wi, 2 = 0. For zero k x or k yi we obtain W r 1 ( 2 ) = 0. 
This agrees with the existence of a single MF mode with 
momentum 7r along a ribbon, identified above. 

From Figure |T] it is apparent that the FS pockets are 
approximately related through the nesting vector Q = 
(0,7r) or ( 77 ,0), so one can ask about the effect of a charge 
(or spin) density wave (CDW or SDW) on the topological 
properties. In the case of a SDW, we keep s~ as a good 
quantum number. The Chern numbers remain the same 
as before, the spectrum remains gapped and similar to 
that in Figure [2] From the point of view of the winding 
numbers, however, there are quantitative changes in the 
case of a CDW. If, for instance, Q = (0, n) then the 
original BZ folds such that 7 t/2 < k y < 7t/ 2. For the 
topological case, t 2 = —0.08, we obtain W\{k x = 7 r) = 
—2 and W\{k x = 0) = 0. This means that the MF still 
has momentum 7r along a ribbon in the xx direction. The 
winding number W 2 {k y = 0) = 1 because the points 
k y = 0,7T are the same under the BZ folding. The MF 
has zero momentum for a ribbon along Q, because of the 
BZ folding. 

In the case of a SDW, the Hamiltonian in equation © 
acquires a term proportional to ro and no longer enjoys 
the chiral properties discussed. 









4 


71 

0 

71 


71 

0 

71 


VA 


71 

71 

+ 

VA 

+ 

A 

J 

O 

c 

0 

0 


O 

c 


dA 


— 7T 

—71 

N 

r\ 

d 

71 

0 

7 T 


71 

0 

71 


kx k x 


FIG. 1: Normal phase FS pockets for model (121-( 131, with 
t 2 = —0.08 (left, topological) and t 2 = 0.08 (right, non- 
topological). 
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FIG. 2: Spectrum for a ribbon geometry. A = 0.1, 1 2 = 
—0.08. The edge modes containing the MF are displayed in 
color. 


C. Induced pairing 


If pairing in the multiorbital superconductor is induced 
by a proximity effect, it would be desirable to know the 
possible matrices r 7 in equations ^ and Q. Such a 
way of engeneering a p + ip supercondutor has been dis¬ 
cussed for single-band systems® and, more recently, for 
s-wave superconductors coupled to Weyl semi-metald^. 
A recent study 31 of the induced pairing symmetries in a 
topological insulator’s surface has been made, in which 
pseudospin degree of freedom is taken into account. It 
was found that a conventional local s-wave superconduc¬ 
tor can still induce spin triplet p+ip pairing, due to spin- 
momentum locking. In our model above there is no spin- 
momentum locking and induced p + ip from a conven¬ 
tional s-wave superconductor is not expected. However, 
if a Rashba-type spin-orbit term is included in Hamilto¬ 
nian then spin-momentum locking will exist and the 
Chern number in equation © will not be affected as 
long as the energy gap is not closed. 


Suppose the Hamiltonian T~L{ A = 0) of the multiorbital 
normal 2D system, in equations 0) and (12), is coupled 
to a single band superconductor H s through a tunneling 


term T, 

H s = ( e (k) — p) &!<. CT fr k ,cr + A a . a ./(k)&j tiCr 6^ k o . / + H.c., 

k,(7 

_ (19) 
T = Y, ^4kA,"+ H - c -- (20) 

k,j=l,2 


If the pairing term in equation (19) is a local spin singlet, 
the induced pairing in the multiorbital layer contains the 
even parity local spin singlet terms (in < 7270 , 1 ) as well as 
the p + ip spin singlet , <4(k)cr 2 T 2 . The above spin triplet 
p + ip term (13) is also induced. The latter is weaker, 
however, and the resulting system is not topological. 

But if the pairing term in equation (19) is a p + ip spin 
triplet, the induced triplet term (13) is stronger than the 
local spin singlet term (cr 2 7o), and the multiorbital sys¬ 
tem has the topological properties discussed above. An 
additional local spin triplet term (in <7 it 2 ) is also induced. 
Other smaller p + ip spin triplet terms with odd (even) 
parity, in d z (k)<Ti7i(3), are also induced. If one of the 
tunneling amplitudes is larger than the other (V\ V 2 ) 

the induced spin singlet term is strongly supressed. 


III. THE ANDREEV PROBLEM 

The Chern number predicts the existence of Majorana 
edge states in the superconductor. Because a Majorana 
operator must be its own antiparticle, it must have mo¬ 
mentum either 0 or it along the edge. In an Andreev 
scattering problem, the Majorana state should appear as 
an Andreev bound state with transverse momentum (ie. 
along the interface) 0 or n. 

In the framework of single band BTK theory, each 
Fermi pocket should contain one zero energy Andreev 
bound state. This is because the pairing function is 
odd, so that the electrons feel a sign change in the gap 
function upon specular reflection at the superconductor’s 
surfac^ 32 } Thus 3 or 4 ABS’s or MF’s are predicted by 
BTK theory, at conflict with the topological properties. 
The single MF observed when only a single FS pocket 
is crossed implies that somehow the two Andreev bound 
states predicted by BTK theory for the pockets (0, 0) 
and (n, 0) should interfere destructively when the inci¬ 
dent electron has transverse momentum k y = 0. 

We consider a N/S boundary along yy axis. In the 
framework of ()\\ r lIAlISl, the incident electron from a sin¬ 
gle band normal metal will split into the two pseudo-spin 
channels of the superconductor, as figure [3] explains. In 
the N side, x < 0, the wavefunction for electrons near 
the Fermi level is exp[ik y y]i/)N where 

ip N (x < 0) = ( J ) jv+ x +b ( l ) e~ ip+x +a ( J ) e*-* . 

( 21 ) 

The momenta p± are close to the Fermi momentum pp 
and are fixed by the energy, E. The amplitudes for elec- 
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FIG. 3: Tight-binding showing the splitting of the incoming 
electron into two pseudo-spin channels of the superconductor, 
similar to a waveguide. 


FIG. 4: Differential conductance at normal incidence (left) 
and for k y = n (right). Model ( |12| ) -(131 with A = 0.1 and 
t 2 = 0.08 (non-topological). Dashed line: normal state differ¬ 
ential conductance. 


tron reflection, b, and Andreev hole reflection, a, al¬ 
low us to obtain the differential conductance as g s = 
1 + |a| 2 — |fe| 2 , whereas the normal state conductance is 
just g n = 1 — \b\ 2 . We shall consider here only k y = 0 
or 7r. The transmitted waves into the superconductor 
are superpositions of wavevectors k , q 31 from two FS 
pockets (see figure |T]) and the wavefunction for x > 0 is 
exp [ik y y]i/js where 

Mx> 0 ) = C^ + e lk+x + D^-e~ lk ~ x 

+ E^+e l!1+x +F^-e- iq ~ x , ( 22 ) 


where each denotes a four-dimensional column eigen¬ 
vector of the BdG matrix, H. The x-components of the 
momenta, fc , , are chosen so that the group veloc¬ 

ity is positive for energy E above the gap. For subgap 
energies, the momenta have a positive imaginary part. 

In QW T t he matching conditions for wave functions 
(21) and (22) at x = 0 are written as 17 18 


^jv(O) <g) 


= V’s(O), 


dkxH N ip N (0) = 


110 0 
0 0 11 


(23) 

dk^Hs^s^x = 0 ). 

(24) 


Here, Hjy denotes a BdG Hamiltonian matrix for the 
normal metal. Interface disorder can be accounted fo]P 
by making the replacement: 1 — b —> 1 — b — 2iZ(l + 
b )pf/p+ and a —> a( 1 — 2iZpp/p + ) in the right-hand 
side of equation (24), and where Z denotes the BTK 
parameter^. 


—0.08) or in the non-topological case (£2 = +0.08). For 
transverse momenta k y = 7 r, equation (25) is verified only 
in the topological case. This means that the quantum in¬ 
terference effects from the two FS pockets effectively ani- 
hilate the two ABS’s that would be predicted by single 
pocket BTK theory. 

One might be tempted to read condition (251 as the re¬ 
quirement that the linear system in equation (23) be ho¬ 
mogeneous and the wave function ips be made to vanish 33 
at the N/S boundary. This would be incorrect, however, 
as it would imply the conductance g s to vanish. On the 
contrary, the conductance g s is finite and independent 
of the dis ord er param eter Z a t the energy value where 
equation (25) is obeyeeP^ |23 | 32 l 


For the non-topological case, the differential conduc¬ 
tance at normal incidence, k y = 0, is shown in figure [4] 
(left) as a function of energy (which in an actual exper¬ 
iment is obtained from the voltage bias). It is seen that 
the quantum interference supresses quasi-particle trans¬ 
mission, g s , as E —> 0 and the effect is even more pro¬ 
nounced as disorder increases. A similar result is ob¬ 
tained at k y = 7T. Features such as peaks and dips are 
visible when the energy E crosses the superconducting 
gaps on the Fermi pockets. 

For the topological case, the differential conductance at 
normal incidence, k y = 0, is shown in figure [5j A similar 
destructive interference is observed at low energy. 

For transverse momentum k y = it , the ABS (MF) 
leaves its imprint on the conductance, as figure [6] shows. 
The differential conductance attains the maximum value 
g s = 2 at E = 0, independent of Z. 


IV. RESULTS 


V. SUMMARY AND CONCLUSIONS 


According to QWT^, the condition for the existence 
of ABS is obtained from the 4x4 matrix A composed of 
the four column vectors 0k+, , </> q +, </> q -, in equation 

(22). The condition then reads 


Det A = 0 . 


(25) 


We checked that condition (25) does not hold for trans¬ 


verse momentum k y = 0 either in the topological (t 2 = 


We have studied a model of a multiorbital 2D super¬ 
conductor which is realized by introducing p + ip Cooper 
pairing, diagonal in the pseudospin channel, in an oth¬ 
erwise TI. Under these conditions, the topological prop¬ 
erties of the TI were shown to become irrelevant. The 
topological properties of the multiorbital superconductor 
depend on the FS in the normal phase. If the FS contains 
several sheets, we have shown that when the number of 
Fermi pockets is odd, the MF is in the single pocket that 
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is traversed by the quasi-particles, while the other pair 
of pockets interfere destructively. This result from QWT 
reconciles the number of MFs with the Chern number, in 
contrast to BTK theory. In addition to the destruction 
of the MFs, the waveguide interference effects also pro¬ 
duce a vanishing conductance at E = 0, when a pair of 
FS pockets is traversed by the quasi-particles. 


FIG. 5: Differential conductance at normal incidence, for 
the same model as in Figure [4] with t 2 = —0.08 (topological). 
Dashed line: normal state differential conductance. 



FIG. 6: Differential conductance for the topological case, 
showing the signiture of the MF. 
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